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I. INTRODUCTION 



The neutrino data show that neutrinos have tiny masses as compared to the electroweak 
scale. This is clear evidence for physics beyond the standard model (SM). The data also 
indicate that the structure of flavor mixing for neutrinos is largely different from that for 
charged leptons. These facts would suggest that, while charged leptons have Dirac type 
masses, the neutrino masses are of the Majorana type. The tiny Majorana masses of left- 
handed neutrinos are generated from the dimension-five effective operators 

£ = |pWV°, (i) 

where A represents a mass scale, are dimensionless coefficients, and <f)° is the Higgs 
boson. After electroweak symmetry breaking, the mass matrix M l J for left-handed neutrinos 
appears as M % J = Qj(0°) 2 /A. As the vacuum expectation value (VEV) (0°) of the Higgs 
boson is 0(100) GeV, the observed tiny neutrino masses (M l J < 0.1 eV) are realized when 
(cy/A) ~ O(10~ 1A ) GeV -1 . It has been an interesting problem how we can naturally explain 
such a small number with less fine tuning. 

If the operators in Eq. ([I]) appear at the tree level in the low energy effective theory, 
A has to be as large as C(10 8 ) - C(10 14 ) GeV for cy being C(10~ 6 ) - 0(1) to describe 
the data. For example, in the tree-level seesaw scenario where right-handed neutrinos are 

n 

introduced, their Majorana masses have to be set much higher than the electroweak scale [![, 
corresponding to the scale A in Eq. (TjQ). Although the scenario is simple, it requires another 
hierarchy between the mass of right-handed neutrinos and the electroweak scale, and in 
addition, physics at such a large mass scale is difficult to be tested at collider experiments. 

Quantum generation of neutrino masses is an alternative way to obtain (cy/A) ~ 
O(10~ w ) GeV -1 . Due to the loop suppression factor, A in these models can be lower 
as compared to that in the tree-level seesaw models. Consequently, the tiny neutrino masses 
would be explained in a natural way by the TeV-scale dynamics without introducing very 
high mass scales. The original model of this line was proposed by Zee where neutrino 
masses were generated at the one-loop level. Some variations were considered UQ, for 
example, by Zee and Babu [3], Kraus-Nasri-Trodden Ma Jo|, and the model in Ref. js]. 
The last three models contain dark matter (DM) candidates with the odd quantum number 
under the discrete Zi symmetry. It must be a charming point in these TeV-scale radia- 



2 



tive seesaw models that they are directly testable at the collider experiments such as Large 
Hadron Collider (LHC) and the International Linear Collider (ILC). 

A general feature in radiative seesaw models is an extended Higgs sector, whose detail is 
strongly model dependent. The discovery of these extra Higgs bosons and detailed measure- 
ments of their properties at current and future collider experiments can give partial evidence 
for the radiative seesaw models. In the literature js-14], phenomenology of these radiative 
seesaw models has already been studied extensively. Such previous works mainly discuss 
constraints on the flavor structure from the current data for such as neutrino physics and 
DM, and also study collider phenomenology of the Higgs sectors 
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Another common feature in radiative seesaw models is the Majorana nature. In order to 
induce tiny Majorana masses for left-handed neutrinos, we need to introduce its origin such 
as lepton number violating interactions in the scalar sector {2, 3| or right-handed neutrinos 
with TeV-scale Majorana masses 4|t6(- When the future data would indicate an extended 
Higgs sector predicted by a specific radiative seesaw model, the direct detection of the 
Majorana property at collider experiments should be a fatal probe to identify the model. 

In this Letter, we study the phenomenology in TeV-scale radiative seesaw models, in 
particular, a possibility of detecting the Majorana nature at collider experiments. We mainly 
discuss three typical radiative seesaw models as reference models; the model by Zee and Babu 



where neutrino masses are generated at the two- loop level [3j, that by Ma with one-loop 
neutrino mass generation and that in Ref. js] where neutrino masses are generated at 
the three-loop level. Typical parameter regions where the data can be satisfied have been 
already studied in each model in the literature. We here study collider phenomenology in 
such typical parameter regions in each model, and discuss the discrimination of these models 
by measuring the details of the Higgs sector and the Majorana nature at the LHC and the 
ILC. 



II. RADIATIVE SEESAW MODELS 



A. The Zee-Babu model 



f. |a 



In the model proposed in Ref. |3|] (we refer to as the Zee-Babu model), in addition to 



singly- charged singlet scalar bosons a; 1 * 1 , doubly-charged singlet scalar fields fc ±± are intro 
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FIG. 1: Feynman diagrams for neutrino masses in the model by Zee-Babu [3] (left), that by Ma [f| 
(center) and that in Ref. [6(] (right). 

duced, both of which carry the lepton number of two-unit, and their interactions are given 
by 



Ant = UL: L Li L )e tJ u + + g' ab (t aR t bR )k ++ ~ V k ++ u-u- + H.c, (2) 

where L L is the left-handed lepton doublet and £ R is the right-handed lepton singlet. The 
matrices and g' ab are respectively an anti-symmetric and a symmetric couplings and the 
lepton number is violated by the interaction with the parameter \x. 

The neutrino mass matrix is generated at the two-loop level via the diagram in Fig. [T] 
(left); 

3 / x \ 2 4^ 

M ii = ^2\T^) Z^fik(ye k 9keye e )fejV 2 h(ml/ml), (3) 
k,i=l \ / w 

where yi [= y/2rrii/v (i = e, /i, r)] are the SM Yukawa coupling constants of charged leptons 
with the masses m 8 and the VEV v (~ 246 GeV), g^ are defined as g„ = g' u and g^ = 2g'^ 
(i 7^ j), m w and are masses of w ± and /c ±± , and 

I 1 ( r ) = -f dx [ dy - l n y ^ ~ y \ (4) 

Jo Jo x+{r-l)y + y 2 x + ry 

where I-i(r) takes the value of around 3 - 0.2 for 10~ 2 < r < 10 2 . The universal scale 
of neutrino masses is determined by the two-loop suppression factor l/(167r 2 ) 2 and the 
lepton number violating parameter fi. The charged lepton Yukawa coupling constants y^ 
(y e <C y^ <C y T < 10~ 2 ) give an additional suppression factor. Thus, any of or g^ can be 
of 0(1) when m w and rrik are at the TeV scale. The flavor structure of the mass matrix is 
determined by the combination of the coupling constants fij and yig^yj. 



The flavor off-diagonal coupling constants and gij induce lepton flavor violation (LFV). 
From the results of /i -> ery, r -> e7 and r ->■ /ry, |/ Mr / re |, 1/^/^1 and |/ Te / e/x | are 
respectively constrained as a function of m u . The data of rare decays of fi — > eee, r — > \i\ie 
and r — > fiee are also used to constrain the combinations \g^ e 9ee\, |flVegw» + fi^fi^e! and 
\9Te9tJ,e + 9riJ.9ee\, respectively, depending on m k . The g — 2 data can also be used to constrain 
a combination of these coupling constants with m u and m^ 1 . 

In the scenario with hierarchical neutrino masses, satisfy f efl ~ f eT ~ The 
typical relative magnitudes among the coupling constants g^ can be g m : g^ T : g TT ~ 1 : 
iTifi/rrir '■ {fTi^/mr) 2 ■ For g^ ~ 1, the neutrino data and the LFV data give the constraints 
such as mfc > 770 GeV and m w > 160 GeV [10j | . On the other hand, the constraints on 
the couplings and masses are more stringent for the inverted neutrino mass hierarchy. The 
current data then gives m u ~ 825 GeV for g^ ~ 1 [10]. One of the notable things in this 
case is the lower bound on sin 2 2#i3, which is predicted as around 0.002 [9]. 



B. Models with TeV-scale Right-handed Neutrinos with a discrete Z2 symmetry 

Similar to the tree-level seesaw model, tiny masses of left-handed neutrinos would also 
come from Majorana masses M^a of gauge-singlet right-handed neutrinos iV£; in the radia- 



tive seesaw scenario 



6J. One simple way to realize the absence of the tree- level Yukawa 
interaction V 1 l <&Nr is introduction of a discrete Z2 symmetry, with the assignment of the 
odd quantum number to iV^ and the even to the SM particles. To obtain the dimension five 
operator in Eq. ([T|) at the loop level, we need to introduce additional Z 2 -odd scalar fields. 
The lightest of all the Z 2 -odd particles can be a candidate of DM if it is electrically neutral. 
The original model of the radiative seesaw model with such a discrete symmetry is proposed 
by Krauss, Nasri, and Trodden JJ], in which neutrino masses are induced at the three-loop 
level. In the following, we consider two variant models of the Kraus- Nasri- Trodden (KNT) 
model. 



1 If we take g ee = 0, then mk is unbounded from the [i — > eee and r — > lee results (£ = e or /i), so that 
relatively light fc ±± (mk ~ 100 — 200 GeV) are possbile. 
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1. The Ma model 

The model in Ref. j^J, which we here refer to as the Ma model, is the simplest radiative 
seesaw model with right-handed neutrinos N^, in which the discrete Z 2 symmetry is intro- 
duced and its odd quantum number is assigned to N^. The Higgs sector is composed of 
two Higgs doublet fields, one of which (H) is Z 2 odd. As long as the Z 2 symmetry is exact, 
the neutral components of H do not receive VEVs. We have one SM-like Higgs boson h, 
and four physical Z 2 -odd scalar states; £° (CP-even), (CP-odd) and ^ as physical scalar 



states. This Z 2 odd Higgs doublet is sometimes called as the inert Higgs doublet [15] or the 
dark scalar doublet jl7| . The LEP II limits are studied in this model in Ref. jl8 |. 

The neutrino masses are generated at the one loop level via the diagram depicted in 
Fig. [1] (center), in which Z 2 odd particles, £° and N^, are in the loop. The mass matrix is 
calculated as 

mk = - f (*) U + _l_ , n A , ( 5) 

^ ) M n% 1 - r a V 1 - r a J' K J 

where hf are the Yukawa coupling constants of V l L 'EN^, M^a is the Majorana mass of 
the a-th generation right-handed neutrino iVjj, A5 = (m| — m| )/v 2 , r a = m^/M^a with 
rriQ = (m^. +m^)/2, where and are masses of £° and respectively. The universal 
scale for neutrino masses is determined by the one-loop suppression factor l/(167r 2 ), A 5 and 
Mj^a. The flavor structure in is realized by the combination of hfh^ / 'Mfja. Therefore, 
for M^a ~ 0(1) TeV, the combination of the coupling constants would be |As|(/if) 2 ~ 10~ 9 . 

In this model, there are two scenarios with respect to the DM candidate; i.e., the lightest 
right-handed neutrino or the lightest Z 2 -odd neutral field (£° or £?). For both cases, 
there are parameter regions where the neutrino data are adjustable without contradicting 



other phenomenological constraints 13| . In this Letter, we consider the case where the dark 
doublet component £° is the DM candidate 2 . When the mass of the DM is around 50 GeV, 
the typical value of A5 ~ 10~ 2 for the neutrino masses gives the mass difference between 
£° and £P about 10 GeV. 3 The relic abundance of such DM is consistent with the WMAP 
data 1191 . 



2 In Ref. [13(, the scenario where the lightest right-handed neutrino is DM is explored. 

3 In order to avoid constraint from the DM direct search results, it is required that | A5 j > 10~ 6 . 
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2. The AKS model 

In the model in Ref. {f5|, which we here refer to as the AKS model, it is intended that not 
only the tiny neutrino masses and DM but also baryon asymmetry of Universe are explained 
at the TeV scale. In addition to the TeV-scale right-handed neutrinos N R (a = 1,2), the 
Higgs sector is composed of Z2-even two Higgs doublets $j (i — 1,2) and ^2-odd charged 
singlets S^ 1 and a i^-odd neutral real singlet rf. Therefore the physical states in the i^-even 
sector are H (CP-even), A (CP-odd), H ± and h (CP-even). 

The neutrino mass matrix is generated at the three-loop level via the diagram in Fig. [T] 
(right), and is expressed as 

v^/ 1 \ 3 (Vt .hf)(y t ./#)(« tan/9) V 
M ij = 2^{Y^) — 3 -Jf~ a I 2 (m H± ,m s± ,m N .,m v ), (6) 

where m,H± , ms±, m N% and m v are the masses of the doublet originated charged Higgs boson 
H ± , S^, N R and r/ , respectively; h" 1 and kv are the coupling constants of lSf a R e R S + and 
H + S~rf } respectively; tan/3 = and 

T( , -4^ 2 f°° , f B 1 (-u;x,y)-B 1 (-u;0,y) \ 2 f z 2 w 2 \ 

I 2 (x,y,z,w) = — 2/ udu< > — —j ■ — r , (7) 

z A — W z Jo I x A \u + z A u + w z J 



where B\ is the tensor coefficient function in the Passarino-Veltman's formalism [25|]. Al- 
though the Higgs sector is rather complicated to make it possible for the electroweak baryoge- 
nesis scenario, the flavor structure is determined only by the combination of hf and m^™ just 
as in the Ma model. The mass matrix has the three loop factor 1 / (167r 2 ) 3 with additional sup- 
pression factor by yi. They are enough to reproduce the neutrino mass scale. Thus, the elec- 
tron associated coupling constants hi' 2 and the scalar coupling k are of 0(1) for rn]f R ~ 0(1) 
TeV. The Yukawa coupling constants are hierarchical as hl ,2 (~ 0(1)) 3> h 1 ^ 2 3> hi' 2 . 

The parameter sets which satisfy the current data from neutrino oscillation, LFV, relic 
abundances of DM and the condition for strongly first order electroweak phase transition 
are studied in Ref. jf], Q. To reproduce the neutrino data, the mass of H ± should be 100 
- 200 GeV. This is an important prediction of the model. In order to avoid the constraint 
from b — y sj, the Yukawa interaction for the doublet fields takes the form of so-called Type- 

n 

X [20j, 4 where only one of the doublets couples to leptons and the rest does to quarks. The 



Type-X is referred to as Type-IV in Ref. Q and Type-F in Ref. 0. 
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physics of the Type-X two Higgs doublet model (THDM) shows many distinctive features 
from the other type of extended Higgs sectors. For example, H and A decay mainly into 
t + t~ when tan/3 > 3 and sin(/3 — a) ~ 1 [20]. There are basically two DM candidates, 77° 
and iVp. The mass of S is strongly constrained by the current data and the requirement 

fin 

for strongly first order phase transition [6j, [14|. The coupling constant of S + S h is required 
to be of 0(1), whose indirect effect appears in the quantum correction to the hhh coupling 
constant as a large deviation from the SM prediction jl^l l^ . As long as kinematically 
allowed, S ± decays via S ± -> #±7/° by 100%. 



III. PHENOMENOLOGY IN RADIATIVE SEESAW MODELS AT THE LHC 

The existence of the extra Higgs bosons such as charged scalar bosons, which are a com- 
mon feature of radiative seesaw models, can be tested at the LHC. Details of the properties 
of such extra Higgs bosons are strongly model dependent, so that we can distinguish models 
via detailed measurements of extra Higgs bosons. In addition, as the (SM-like) Higgs boson 



h is expected to be detected, its mass and decay properties are thoroughly measured [29]. 
The radiative seesaw models with a DM candidate can also be indirectly tested via the in- 
visible decay of h as long as its branching ratio is more than about 25% for = 120 GeV 



with C = 30 fb^ 1 
literature 
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he phenomenological analyses at the LHC in each model are in the 
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20N24j] . We here review some remarkable features. 
At the LHC, u^ 1 and fc ±=t in the Zee-Babu model can be produced in pair, via the Drell- 
Yan s-channel processes qq — » u + ui~ and qq — > k ++ k . The direct detection of k ±±: can 
be a signature for this model. The cr(qq — > k ++ k ) is around 0.1 fb for ~ 800 GeV. 
If fc ±± mainly decay into Afi (or 4e), the number of the signal event with £=100 fb -1 [31 1 



is enough for the discovery for nik < 800 GeV (9]. The doubly- charged Higgs bosons are 
however also predicted in the models with complex triplet scalar fields, A = (A ±=t , A ± , A ). 
The gauge coupling W ± A ± A TT induces the single doubly- charged Higgs production q q' — > 



\y±* _v. A±±A=F ; whose cross section is comparable to that of qq — > 7*, Z* — > A ++ A [32]. 
The absence of the gauge coupling W ± u ± k TT is an important distinctive feature of the Zee- 
Babu model with gauge singlet doubly-charged Higgs bosons from the triplet model. The 
singly-charged Higgs boson ^ would be more difficult to see the signal at the LHC because 
the final state from u + u>~ is £ + £~ plus a missing energy. 
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In the Ma model, if the mass of the DM candidate satisfies m 



16 



[ri < m h /2, the SM Higgs 



17j. For m (r = 50 GeV and 



boson with rrih < 2mw decays mainly into the DM pairs 

= 60 GeV, the branching ratio of the invisible decay h — > reaches to 70% around 
rrih ~ 120 GeV, which can be observed at the LHC. The h — > 77 mode is suppressed by 



an additional contributions from [17]. In Ref. [17], the discovery potential of the dark 
scalar doublet is also analyzed by pp — >• for the benchmark points, m^ r = 50 GeV, 
= 60 - 80 GeV, and m^± = 170 GeV. The decay branching ratio of the CP-odd dark scalar, 
# -»■ Z*^ -»■ #£ r , is about 0.09 (0.07) for - m & = 10 (30) GeV. A signature ii plus a 
missing energy in the benchmark scenario may be discovered by the optimal cuts. 

In the AKS model, the invisible decay of the SM-like Higgs boson h — > rj rj can also 
open if kinematically allowed. For the typical scenario in Ref. [0, Q], the branching ratio of 
the invisible decay can amount to B(h — > rfrf) ~ 36 (34) % for m v = 48 GeV, rrih = 120 
GeV and tan/3 = 3 (10), so it would be testable at the LHC. As the coupling of hS + S~ 
is strong, the partial width of Y{h — > 77) deviates from the SM prediction. The lepton 
specific decays of extra Higgs bosons A, H and H ± are discriminative feature of the Type-X 
THDM [2oJ-|24|. The dominant decay mode of H (A) is H t + t~; B(H(A) -»■ t + t~) ~ 1 
for tan/3 > 3. The decay into is suppressed by a factor of {m^/m T ) 2 . These neutral 

bosons can be seen by gg h, A, H ^ t + t~ (/i + /i~) |20|, [22j, [24j. The doublet originated 



charged Higgs boson H^ 1 (as well as extra neutral ones) is as light as 100 - 200 GeV, so t 



hat 



the property of the Type-X 



20 



23] 



TIDM can also be tested by pp — > AH ± — » t + t^t ± i i t 
and pp — > AH — > 4r 20j, [2JJ] . On the other hand, the mass of the Z2-odd charged Higgs 
bosons S ± is around 400 GeV in the typical scenario in Ref. Q. They are produced 
in pair via the Drell-Yan process, and decay as S + S~ — > H + H~rfrf — > T + T~vi>rfrf . The 
event rate is about 0.5 fb for ms± = 400 GeV when tan/3 > 2. Separation of the S + S~ 
signal from the H + H~ event and also the SM backgrounds seems to be challenging. 

At the LHC, via the physics of extra scalar bosons such as (singly and/or doubly) charged 
Higgs bosons and CP-even Higgs bosons, the structure of the extended Higgs sector can be 
clarified to some extent. In addition, the invisible decay of the SM-like Higgs boson and the 
mass spectrum of the extra Higgs bosons would give important indication for a possibility 
to a radiative seesaw scenario. However, although they would be a strong indication of 
radiative seesaw models, one cannot conclude that such Higgs sector is of the radiative 
seesaw models. In order to further explore the possibility to such models, we have to explore 
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the other common feature of radiative seesaw models, such as the Majorana nature. In the 
next section, we discuss a possibility of testing the Majorana nature at ILC experiments. 



IV. PHENOMENOLOGY IN RADIATIVE SEESAW MODELS AT THE ILC 

At the ILC, properties of the Higgs sector can be measured with much better accuracy 
than at the LHC, so that we would be able to reconstruct the Higgs potential in any extended 
Higgs sector if kinematically accessible. Invisible decays of the Higgs boson can also be tested 
when the branching ratio B(h — > invisible) is larger than a few % [331 • Furthermore, the 
Majorana nature in radiative seesaw models; i.e., the existence of TeV scale right-handed 
Majorana neutrinos or that of lepton number violating interaction, would also be tested at 
the ILC. 



A. Electron-positron collisions 

In the pair production of charged scalar bosons at the e + e~ collision, which appear in 
the radiative seesaw models (u + u~ in the Zee-Babu model, in the Ma model, and 

S + S~ (and H + H~) in the AKS model), there are diagrams of the t-channel exchange of 
left-handed neutrinos or right-handed neutrinos in addition to the usual Drell-Yan type s- 
channel diagrams. The contribution of these t-channel diagrams is one of the discriminative 
features of radiative seesaw models, and no such contribution enters into the other extended 
Higgs models such as the THDM. 5 These t-channel effects show specific dependences on 
the center-of-mass energy \fs in proportion to logs in the production cross section, and 
enhances the production rates of the signal events for higher values of y/s. The final states 
of produced charged scalar boson pairs are quite model dependent but with missing energies; 

e + e _ ->■ oj + uj~ ->■ , [The Zee-Babu model] (8) 

e + e - ^ ^ W^W-^e^r .I.IJ.l[.ljlii!L)CC , [The Ma model] (9) 
e+e~ -> S + S- -> H+H-rfrf -> t%t^ V l v l ^ , [The AKS model] (10) 

5 In the minimal supersymmetric standard model (MSSM) selectron pair production can have similar t- 
channcl contributions (the Bino exchange). In such a case, the final state would be something like a e + e~ 
pair plus a missing energy. Therefore, we can discriminate it from the radiative seesaw models. 
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FIG. 2: The differential cross section of e + e~ — > uj + lj~ — > /i~T + (+missing energy) as a function 
of the angle between the outgoing muon and the beam axis in the Zee-Babu model for y/s = 1 
TeV. The rate of /i~r + (+missing energy) from main background e + e~ — > W + W~ is also shown. 

where underlined parts in the final states are observed as missing energies. 

The Zee-Babu model 

In the Zee-Babu model, the decay branching pattern of uj^ is determined by the relative 
magnitudes of the coupling constants /y. As a reference scenario, we take a parameter set 

mu] = 300 GeV, m k = 1200 GeV, /i = 800 GeV, 
U = fer = 0.013, / MT = 0.027, 

g ee = 0.17, = 1.8, g TT = 0.0061, g e/1 = 5.7 x 10~ 5 , g er = 0.011, 9flT = -0.081, (11) 

which satisfies the neutrino data for the normal mass hierarchy. In this scenario, the rate 
in final states £ + £~ = e + e~, e ± fi Zf , e ± r T , /U ± r =F , r + r _ is given by 2, 13, 13, 19, 36, 

19, respectively. In Fig. 2, the differential cross section da /d cos 9^ for the signal e + e~ — > 
u + u~ — > n~r + (+missing energy) is shown for y/s = 1 TeV as a function of cos 6^, where 9^ 
is the angle between the outgoing muon and the beam axis. The main background comes 
from the W pair production, which is also plotted. The angle cut (e.g. cos^ < —0.5) 
improves the ratio of the signal and the background. 6 

6 Although in this Letter we mainly discuss the case where mk is at the TeV scale, we just comment on the 
case of lighter k . In such case, the pair production of k ++ k can be a clear signature of this model, 
whose signal is the like-sign dilepton pairs with opposite direction 0, Q . 
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FIG. 3: The jets invariant mass distributions of the production rates of the signal in the Ma 
model at ^/s = 500 GeV. left : The di-jet invariant mass M(jj) distribution of the signal 
e + e~ — > £ ' £ — > jj/J>v£,r£r f° r m £± = GeV. right : M(jjjj) distribution of e + e~ — > — > 
— > jjjjCr f° r m £± = 150 GeV. In addition to the rate from the signal process, those 
for main backgrounds are also shown. 

At current and future LFV experiments, the coupling constants fij and can be further 
tested via the LFV rare decays such as £ — > £'7 and £ — > £'£'£". The same operators as in 
£~ — > £' T e~e ± would also be tested directly at the ILC via e ± e~ — > £~£! ± . In the scenario 
in Eq. ( 1TT1) . we estimate that a(e + e~ — > /i ± r T ) ~ 5 fb for i/i = 1 TeV. When we take 
g ee = 0.4, g er = 0.01, ink — u — 1-2 TeV and m w = 400 GeV, which also satisfy all the 
current data, we obtain a(e + e~ — > r ± e T ) ~ 0.76 (1.7) fb for y/s = 500 GeV (1 TeV). 

The Ma model 

In the Ma model, the coupling constants (a = 1,2) 7 are strongly constrained from 
neutrino data and LFV data. As a typical choice of parameters, we consider 8 

mg r = 50 GeV, m& = 60 GeV, m^± ~ 100 GeV, = 171^ = 3 TeV, 

A 5 = -1.8 x 10~ 2 , h%, h°, K ~ !0" 5 , (12) 

in which the normal neutrino mass hierarchy is realized. Because hi are very small for a 

7 Here we consider the minimal case of two generations for the right-handed neutrino. 

8 The relatively large mass difference between ^ and £° i implies a significant deviation from the custodial 
symmetry in the Higgs sector, which affects the allowed mass of the SM like Higgs boson h. The 
larger m,h is favored for larger mass difference of m^± — m^ r i . 
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TeV scale raw™, the contribution of the t-channel diagrams to the signal e + e~ — > is 
much smaller than that from Drell-Yan type diagrams. For most of the possible values of 

and m N a which satisfy the LFV and the neutrino data, the contribution of the t-channel 
diagrams is negligible. The production cross section of a charged Higgs pair is therefore 
similar to that in the usual THDM: about 92 (10) fb for m 5 ± = 100 (150) GeV at y/s = 500 
GeV. The produced decay into W ±w C°r 9 

In Fig. 3 {left), we show the invariant mass distribution of the di-jet jj of the production 
cross section of the signal, e + e~ ->■ ->■ -»■ jjfiv£r& f° r m £± = 100 GeV. 

The main backgrounds come from WW. The jjfifi events from ZZ, 77, and Z7 can also 
be the backgrounds. A factor of 0.1 is multiplied to the rate of the jjfifi backgrounds for 
the miss-identification probability of a muon. The signal is significant around M(Jj) ~30 
GeV. The invariant mass cut (such as 15 GeV< M(Jj) < 40 GeV) is effective to reduce the 



backgrounds. For the numerical evaluation, we have used a package CalcHEP 2.5.4 34]. 

For m^± > mw + rn^ r , on the other hand, the signal W + W~^^ can be measured by 
detecting the events of four jets with a missing energy. The main background comes from 
W + W~vv and it. By the invariant mass cuts of two-jet pairs at the W boson mass, the 
biggest background from WW can be eliminated. In Fig. 3 (right), we show the invariant 
mass distribution of jjjj of the production cross sections of the signal and the backgrounds 
without any cut. A factor of 0.1 is multiplied to the rate of tt background, by which the 
probability of the lepton from a W that escapes from detection is approximately taken into 
account. The signal is already significant. The invariant mass cut {M(jjjj) < 300 GeV) 
gives an improvement for the signal/background ratio. 

The AKS model 

For the AKS model, we take a typical successful scenario for the neutrino data with the 
the normal mass hierarchy, the LFV data and the DM data as well as the condition for 
strongly first order phase transition Jfj[ Jjj]; 

m v = 50 GeV, m H ± = 100 GeV, m s ± = 400 GeV, m N i = m N 2 = 3 TeV, 

Ft Ft 

h 1 e = hl = 2^>hl,hl^>h 1 T ,h 2 T , k~0(1), sin(/3-a) = l, tan/3 = 10. (13) 



9 The production can also be interesting. The final state should be two jets (or dilepton) plus a 
missing energy. The cross section for e + e~ — > — > £r£rii is about 40 fb at y/s = 500 GeV. 
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FIG. 4: left : The cross sections of the signal, e + e~ — > S + S~ — > t + t~ (+ missing energy), in the 
AKS model as a function of the collision energy y/s. right : The differential cross section of the 
signal for -^/s = 1 TeV as a function of the angle of the direction of the outgoing t~ and the beam 
axis of incident electrons. In addition to the rate from the signal, those from backgrounds such as 
t + t~, t + t~vv and H + H~ are also shown. 

Because hi' 2 ~ 0(1), the contribution from the t-channel exchange diagrams to the 
production cross section of S + S~ dominate that from the Drell-Yan diagrams The 
cross section is about 87 fb for m s ± = 400 GeV at y/s = 1 TeV. As the decay branching 
ratio of — > H rf is 100% and that of H^ 1 — > t^v is also almost 100% because of the 
Type-X THDM interaction for tan/3 = 10, the final state of the signal is T + T~Vvrfrp with 
almost the same rate as the parent S + S~ production. The main SM backgrounds are t + t~ 
and t + t~Vv. The pair production of the doublet like charged Higgs boson H + H~ can also 
be the background. As the signal rate dominantly comes from the t-channel diagram, it 
becomes larger for larger ^fs, while the main backgrounds except for rrv e v e are smaller 
because they are dominantly s-channel processes (Fig. 4 (left)). At = 1 TeV, the rate of 
the signal without cut is already large enough as compared to those of the backgrounds. It is 
expected that making appropriate kinematic cuts will improve the signal background ratio 
to a considerable extent. The \fs scan will help us to confirm that the signal rate comes 
from the t-channel diagrams. Fig. 4 (right) shows the differential cross section of the signal 
at ^/i = 1 TeV as a function of cos 9 T - , where 9 T - is the angle between the direction of the 
outgoing r~ and the beam axis of incident electrons. The distribution of the background 
from tt is asymmetric, so that the angle cut for larger cos 9 T - reduces the backgrounds. 
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B. Electron-electron collisions 



As already stated, the ILC has a further advantage to test radiative seesaw models via 
the experiment at the e~e~ collision option, where dimension five operator of e~e~0 + + , 
which is the sub-diagram of the loop diagrams for neutrino mass matrix. This direct test of 
the dimension five operator is essential to identify the radiative seesaw models. 

The Majorana nature in the Zee-Babu model is in the lepton number violating coupling 
constant fi of k ++ u~u~, which generates the dimension five operator of e~e~u + u + at the 
tree level via the s-channel k exchange diagram. The cross section of e~e~ — > u~u~ is 
given by 



->----) = ^V 1 4^ ^n - (14) 

where the total width Y k of /c ±=t is computed as about 168 GeV in our scenario in Eq. ljTTI) . 
On the other hand, in the Ma model and the AKS model, the operator comes from the 
t-channel right-handed neutrino exchange diagram. The cross section is evaluated as 

2 



o[e e — Y (ft (p ) = / dt 



Cmax j 



1287TS 



n 



(IclW) ^ + 



,2 



/ ,11° I ""H%J ~ 2 _ 

\ L Tfl \JOL Vi Tfl JsJOl 



(15) 



where n is the number of generation of right-handed neutrinos, (f)~ represents the Z 2 -odd 
charged scalar boson £~ in the Ma model and S~ in the AKS model. The constants c a 
represent or in the Ma model or the AKS model, respectively. We note that due to 
the Majorana nature of the t-channel diagram, we obtain much larger cross section in the 
e~e~ collision than in the e + e~ collision in each model assuming the same collision energy. 

The mass matrix of left-handed neutrinos is generated at the one, two and three loop 
levels in the Ma model, the Zee-Babu model and the AKS model, respectively. Therefore, 
the coupling constants can be basically hierarchical among the models, so are the cross 
sections. For the typical scenarios in these models, the cross sections are shown in Fig 5. 
The rate of uj~uj~ in the Zee-Babu model can be larger than several times 100 fb for 800 
GeV < y/s < 1.5 TeV. It becomes maximal (several times pb) at y/s ~ m^, and above that 
asymptotically reduces by l/s. The maximal value of the cross section is sensitive to the 
value of g ee and //. In the parameter sets where these coupling constants are smaller the 
cross section becomes smaller. The signal should be like-sign dilepton with a missing energy. 
On the other hand, in the Ma model, production cross sections of e~e~ — > are smaller 
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FIG. 5: The cross sections of like-sign charged Higgs pair productions in the Zee-Babu model 
{bj~uo~) and in the AKS model (S~S~) are shown as a function of the collision energy i/s. The 
parameters in the Zee-Babu and the AKS model are taken as in Eq. (|lip and Eq. (|13p . respectively. 



than 10~ 4 fb because the coupling constants hf are very small in the parameters in Eq. (|T2|) . 
Allowing some fine tuning, hf may be at most 0.01 for heavier N^. In any case, the cross 
section of e~ er — > is smaller than 10~ 3 fb. Hence, most of the successful scenarios in 
the Ma model the process e~e~ — > is difficult to be seen. In the AKS model, the cross 
section of e~e~ — > S~S~ is large, and its value amounts to about 15 pb at y/s = 1 TeV in 
the scenario given in Eq. ( !T3|) . Above the threshold, the magnitude of the cross sections are 
not sensitive to a/s, so that even if m s± would be at the TeV scale, we might be able to test 
it at future multi-TeV linear colliders, such as the Compact Linear Collider 35]. Because 
i?(S' ± — > T]°H ) ~ B(H ± — > r^v) ~ 100 %, the signal should be T~T~vvriri with almost the 
same rate as long as m s± < m^a. 

The background mainly comes from W~W~v e v ei and the cross section is about 2.3 fb (22 
fb) for -y/i = 500 GeV (1 TeV). The branching ratio for the leptonic decay of W bosons is 
30%, so that the rate of the final state iPvvvv is at most 2 fb or less. Therefore, the signal 
in the AKS model and in the Zee-Babu model can be seen. 

Apart from the TeV-scale radiative seesaw models, there are many models with lepton 
number violating interactions or right-handed Majorana neutrinos. Atwood et al. have 
discussed the signature of heavy Majorana neutrinos in the model without Z 2 symmetry via 
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charged Higgs pair production at e + e~ and e~e~ collisions [36j. 10 In supersymmetric models, 
Majorana particles also appear, and their effects also give similar t-channel contributions to 
the above models in the slepton pair production through the gauge couplings, e~e~ — > e~e~, 
whose cross section is of (9(100) fb. The final state would be e~e~x°X° f° r example. 

The e~e~ collision experiment is useful to test the Majorana nature of radiative seesaw 
models such as the Zee-Babu model and the AKS model via like-sign pair production of 
charged scalar bosons. The cross section can be significant and hierarchical among these 
models. The signal can be observed as a model dependent final state, by which we can 
discriminate the models. Although in this Letter we did not explicitly discuss the KNT 
model that also contains TeV-scale right-handed neutrinos N^, we found that the cross 
section of e~e~ — > S 2 S2 is the Z 2 -odd isosinglet charged scalar boson) is very small 



because the coupling constants for 6rN^S 2 (a = 1,2) are tiny [12 ]. 



V. CONCLUSION 

We have discussed general features of TeV-scale radiative seesaw models. They are char- 
acterized by an extended scalar (Higgs) sector and the Majorana nature. We have mainly 
discussed the concrete models with neutrino mass generation at one- loop (the Ma model), 
two-loop (the Zee-Babu model), and at three-loop (the model in Ref. [6[). Various phe- 
nomenological aspects of these models have been discussed especially in experiments at the 
LHC and at the ILC We have found that, while the extended Higgs sector can be explored 
at the LHC, the Majorana nature of the models can directly be tested at the ILC via the 
pair production of the charged scalar bosons at the electron-positron and electron-electron 
collision experiments. The detailed realistic simulation has to be done elsewhere. 

Acknowledgments 

We like to thank Ernest Ma for useful discussions. The work of SK was supported, in 
part, by Grant-in-Aid for scientific research (C), Japan Society for the Promotion of Science 



10 



Recently, charged Higgs pair production at e e collider is studied in a specific model with three TeV-scale 



right-handed neutrinos and four Higgs doublets in Ref. 37(. 
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